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ABSTRACT 


A two-dimensional omega equation is derived by 
combination of the vorticity and thermodynamic equations, 
The desired omega is then taken to be the logarithmic average 
in the 1000-700 mb layer. A diabatic term, after Laevastu, 
for oceanic areas only is included to deduce the empirical 
temperature and vapor-pressure changes associated with 
sensible and latent heating in the maritime layers. Over 
both continental and oceanic areas a frictional vorticity 
sink is included in order that excessive energy cannot be 
generated over the ocean. Among other novel features is the 
use of the Holl static-stability parameter which affords 
vertical consistency wit analyses prepared by Fleet 


Numerical Weather Facility. 
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1. Introduction 

Little practical use has been made of vertical motions, 
We dp/dt, in short-range forecasting due to the complicated 
and lengthy computations required. Vertical motion computa- 
tions are usually the by-product of multilevel baroclinic 
models and single values must be extracted only after going 
through the entire three-dimensional procedure, In this 
paper a two-dimensional vertical motion equation is developed 
using vertically integrated parameters with an assumed vertical 
motion profile. The layer-mean omega is especially useful 
for short-range thickness forecasts and the 1000-700 mb layer 
has been selected for study. Derivation of the(J -equation 
is similar to that by Thompson [18| except a diabatic heating 
term providing a mechanism for development over oceanic 
areas has been included. Due to the empiricisms employed 
the non-elliptic conditions mentioned by Pedersen [16] до 


not arise. 


2. The Basic Omega Equation 

The diagnostic omega equation is derived by standard 
means from a vorticity equation and a thermodynamic equation 
which retains the diabatic term, Equation (1) 18 the time 
derivative of the first law of thermodynamics in (x, y, p, t) 


coordinates, 


ет. МуТт , 1288 = о (1) 
ot Ө др С» 


Substituting, 7 --195. , which comes from the hydrostatic 
Р 
assumption and the ideal gas law; then dividing by, -g/R, 


yields equation (2). 


AF) Pa) Mo Сана ча M (2) 


Operating on equation (2) with the Laplacian gives the final 


form of the thermodynamic equation as shown by equation (3). 


ray әв) “(уе >s) | Vou. -AVQ (3) 

bp ХАС + (5%) М Ра OCP 

Here, СА « — Rp108, is the Holl stability parameter to be 
ЭЭР 


further discussed in section (li). 


The vorticity equation in pressure coordinates is shown 


by equation (lh). 


95 + М: 3-8), 25-68) = 


(г.р, (дады сад 
Le er T 


As presented by Thompson [18 | ana numerous other writers 

the last term of the left side оѓ (1) is approximately equal 
to that on the right side, and the two terms are henceforth 
deleted, Then, a usoful form of the vorticity equation is 
arrived at (equation 5) by making the geostrophic assumption 


for vorticity and for velocity, $a - 10, Має, апа ру 
taking the logarithmic pressure derivative, 


12 


аж. 2 18,1) - Mo o 
dlp NBE Hp 4 oA» ap 
tere, Y) Se, is the absolute vorticity. Next subtract 


(5) from (3), and the result is the omega equation: 
“(- X5 , | 
уа) + БОЈА m 
-ыУ/Ут/да) ЫҚү 9 
т T) ге 


3. Vertical Distribution of Pressure and Omega 
Instead of the distribution of pressure indexes normally 
considered, that employed here is given by, | 
Th. = Ps (2%) = 27% 
о о 
where Y) is an arbitrary pressure level. The 1000-500 mb 
layer is divided into sub-layers as shown by Figure (1). 
The resultant pressures, while essentially logarithmic, bear 


values nearly identical to those of the mandatory levels. 


Ш ааа Py = Py2 = 500.0 mb 
3 --------- Pa = р,274 = 594.6 mb 
2 --------- Pa = P2 = = 707.1 mb 
1. nest Р, з p 27% = 840.9 mb 
ше: 56 Ра = Бг“ = 1000.0 mb 


Figure 1. Pressure Distribution in the Vertical 
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The values of L0» dP/dt, шаш n y аб == 
sure coordinates, is assumed to vary parabolically in the 
vertical with a profile defined by equation (m which 
extends nearly to the 500 mb level. 


2 
(о = 1025 + В АР | (7) 
Р з 

For boundary conditions it is assumed that O-Oat p-p, 
(the kinematic boundary condition at the assumed surface of 
the earth, p,). Also for convenience it is assumed that at b, 
9%ь-0 , i.e. that the 500-mb divergence is zero. Then, 
upon differentiation of equation (7) with the 500-mb diver- 
gence taken to be zero, equation (8) results. It should be 


noted that if the value ( двубр)= О is assumed zero at p-p, 
the value of B is altered by only 57, 


òw) 2A (aP) 2 Ë B9(L PL o 
(36) | Po op Po Р ЭР Po РАЂА | 


ог, г4( (Ы, аа = ур Ч) С) wo 
Po Гч Ч % м 

Therefore, B=1.3863A, and the omega profile in terms of A is 

given by equation (9). 


en | (ln 15463 (lu FL (9) 
RÀ * Р 


O 
Since, Wo =0, the ratio between CO, and C94, defines the 


omega profile (equation 10) for the 0-2 layer (see Figure 2). 
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(10) 











ЦО —> 
Da = 1,0143 CO | 


Figure 2. The Omega Profile 


For the pressure-range considered (see igure 2), the layer- 
logarithmic-mean, С) =0.950G),, s0 that ә, will be used Гога), 
Note finally, the profile assumed here applies only бо the 
large-scale adiabatic, frictionless component of the vertical 


velocity (that is terrain irregularities are considered 


absent), 


li, Vertical Integration of the Omega Equation 

The Holl stability parameter [11] is used since it is one 
of the vertical-consistency requirements used by FNWF (Fleet 
Numerical Weather Facility) and this study uses FNWF processed 
data. The Holl parameter is a modification of the standard 
stability parameter used here, as shown by equations (11), (12), 
and (13). 
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Он = Кре = де“ ~ = Rp Br - a Эр 


Letting, F] = — A. gives, 
Он = Ка дуг. 4- 12.21 АГ а 
Ср 204% op 


Equation (12) may now be finite differenced over the 0-2 layer. 


Gas 5. R (2.-2.)_ (Ta- -T2) (13) 
lue Cp 
Equation (13) xa ip the Holl stability parameter in 
finite-difference form. In effect, On , by (13) gives the 
logarithmio- pressure Plaverago € over the 052 layer. Тие 
use of os is compatible with the assumed Lorem Г pressure 
distribution and will be used throughout the analysis which 


follows, 


The omega equation (6) will now be integrated term by 
term over the isobaric layer (Po, P2). From (6), one may 


rewrite the QJ -equation as: 


A) 


(0) 
туо) + Fp da = ee e) 
А (E), 
Р геј y Evo 


Whose parts (A),...., (E) will now be discussed individually. 
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(A) The first term is approximated by substituting 


layer mean values; Vs 2 i zi E 29 





(B) With the mathematical relation, рід - 209. - 29), 


MEy- 


whose center-finite-differenced form over the 0-2 layer 


the second term becomes, 


becomes: 


СУ, We- ZwtWo y -(до (121) 
(Elu P) (ҖУДА 

Then, using the specified vertical motion profile the rela- 
tion®, = 1.71430, is substituted. However, it is desired 
to introduce lower boundary effects due to friction and 
terrain by letting (dos COLo , the COLo term then being pro- 
determined and placed on the forcing function side of the ц)- 
equation, Expression (15) below is the final form of (il) 


including effects of the lower boundary. 


(54. Ре, Y" Ju Pe / >, 


LFN бы со _ - 


«УАТАЕРАУЯ 


The lower boundary, Wo, will be discussed in section (5). 


Term (B) = e 023570, + | 2142 ¿O 


(C) Term (C) in the geostrophic-diagnostic model repre- 
sents geostrophic advection of absolute vorticity. Averaging 
in the vertical with respect to the logarithm of pressure leads 


to equation (16): 


ЕГ 


Po Po 
d) Де, п». 24 21(=,0)2 4 a6 
E Pa pif 
The derivation of equation (16) makes use of the mean value 
theorem, so that p = Py is represented by Py and may be re- 
moved from within the integral. Equation (16) then becomes, 


Ра 
Ке? = 22 Р X T (20,0) - ты (е...) (17) 


(D) Term (D) involves the geostrophic advection of 


thickness within the logarthmic pressure integral. Its 
value is ы approximated by equation (18): 


15, Чи у. we Ма 2 је а аа 


the right side of which gives equation E 


ра У JE, Az) = did 23 V “ы, Za- =.) (19) 
fe, ОМ, 
Equation (19) shows that the geostrophic advecting wind may bo 


taken arbitrarily as the level (1) wind. 


(E) The diabatic heating term will be vertically 


integrated in section (6). 


The integrated form of the resulting omega equation, subject 
to the lower and upper boundary conditions (at Pe and рч.) 


becomes, 
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wow ) "p € [2049024 _ 2 RN 


Me PE гра Те, ^ aea) 


ата w 


Equation Lm) is a Helmholtz-type equation in the variable, 
TO, . This ¡ardupias, Ou co, » precludes making the usual 
simplifying approximation of most three-dimensional models, 


See for example, Haltiner et SE 
V (Ç u) ау + «ev'e 4 zVr- co 4 0 Wo 


The final solution of equation (20) is to be divided by Gy 
leaving CO,z (о , where W is the mean vertical motion which 


approximates the 850 mb vertical motion. 


5. The Lower Boundary Condition 

Vertical motion at the lower boundary Wo, is the sum of 
a terrain anà a friction term, Wrei = COT Wr. 
Here, (Or , the terrain-effected vertical motion may be de- 


scribed by equation (21) from Berkofski and Bertoni [2] ° 


= -Ра Ут -VZr (21) 
Terrain height, Zr, is a smoothed field used operationally 


by FNWF, The wind velocity at terrain level, У-- s Will be 


arrived at by an objective procedure to be described below, 


29 


Тие Ргісбіоп бегт(Фр, as shown by equation (22) is a 


simplified version of Cressman's formula used by Haltiner 
et a1| 9| . | 


(Or таны 6% %- (22) 


Сыз» the geostrophic drag coefficient derived by Cressman | 7 |, 
the field of which is in regular use at FNWF. The field of 
С» is a set of constants, one for each (1) grid-point in the 


Northern Hemisphere, 


Since the T subscripts infer application at terrain height, 
values are used that approximate the particular terrain heights 
involved, This is done by dividing the O-2 layer into three 


terrain-height dependent cases, 


masse 3: 2200 m <= =+ 
use P+, values 
Case 2: 900 mgZ,£*2200 п 


use P, values 


Case 1: sfesZr<900 m 


use Pe values 





Figure 3, Three Cases of the Lower Boundary 


Division of the 0-2 layer as shown by figure (3) is 
dependent upon the gradient level above terrain, since Vr is 
approximated by some geostrophic wind throughout. As a basis 
for the selection of the particular geostrophio wind, Wg , note 


that the gradient level in a neutral atmosphere occurs at 


20 


4% 


approximately 700 m above terrain level for a wide range of 
surface roughness, as shown by Blackadar | 3 | . "IhersTore, 
vertical velocity at the lower boundary is case GL), if the 
terrain height is 900mor lower. This values corresponds to 
a "standard atmosphere" height which is 700m ог more below 
the (1) level, and (0) level parameters are then used to 
compute the lower boundary. Case (2) is to be used if the 
terrain height Zr lies within a gradient-level range of p, 
but not Р, . Case (3) occurs for all terrain heights within 


700 m of P+ ° 


Equation (23) is the terrain-induced vertical velocity at 
the lower boundary after application of the geostrophic 


assumption, 


vb. (23) 





г 3 
29) = -fh Рр ТО, гат) -h £, Ó, Ix Dx 


Here, the Y subscript refers to the levels O, 1, and 2. 
Densities are standard atmosphere values for the respective 


levels, 


6. The Diabatic Heating Term 

Diabatic heating in this model results from an exchange 
of heat between the sea surface and the atmosphere, The ex- 
change is divided into three cases: neutral, convective, and 
inversion, the division being dependent upon the sea-air 


temperature difference, Ти - ЦЭС 


` 


The neutral case (0% а 8909 results in no net 


transport of heat across the air-ocean interface and the 


aL 


diabatic heating term ís arbitrarily taken as zero when this 


condition prevails, 


The convective case ( AS B o pesults in a net heating 
of the atmosphere, The convective atmosphere is assumed to 
consist of a dry adiabatic lapse rate to the condensation 
level and a moist adiabatic lapse rate above. This particular 
model of the convective case is discussed by Burke ЕЙ in hís 
paper on the transformation of cP to mP air, Heat gain by the 

atmosphere is manifested by a layer-mean temperature gain 
induced by a one-hour trajectory of surface air over a warmer 
sea-surface, and the layer temperature change is shown by 


figure (4), 


Р. 
Im 


7 = — Fe " 
Fe pas N Ха 
hh NEM 


Po 


eat тчййыйрчугучташё: TO GIA ne ул. Ac A. 


Figure Ц. Model of the Convective-Atmosphere Modification 


Heating is composed of: (1) sensible heat from the sur- 
face to the condensation level and (2) latent heating above 
this level. This result follows by an analysis of the Burke 
model and its implications reguarding heat transport. Martin 
ЕЛ has shown that the large-scale transports of heat have 
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already been included in those terms of equation (6) in 
which Q does not specifically appear. Thus О/бр, for use 


in equation (6) is given by the local temperature change. 


2” са г. (сь АМА + (с Als (21) 
SENSIBLE At /батект 
Equation m is the diabatic heating бегі. (5) vovorve’ 


~ i Mira i | т | а 
ае ата О Са 


„брза „Врх "| Др) – (419) m 
Су ДЕ ДЕ 
Integrating logarithmically in the vertical yields the 


final form of IR diabatic term as shown by equation (26). 


“lp aria Ku [гат ATA Sh, 
go Где 


у" ret 


= — LA У | e- Рата + Pa- Ре) АТА (26) 


In each of euo 2h, 25, and 26, the ME 5 
the finite-difference version of the local time derivitive, 
while, Pa , represents a convective-condensation level 


pre sure. 


According to the model depicted in figure (1), the modi- 
fication is largely determined by the static stability which 
in turn largely depends on Sh . As an assumed infinite source 
of heat, the ocean will first modify a thin surface layer of 

air and then, due to convective activity, will distribute the 
modification (temperature change) throughout the 0-2 layer. 


An overestimation of heat exchange should be expected with 


23 


this particular process since convection is not instantaneous 
This overestimation may be partially balanced by heat lost 

from the top of the 0-2 layer. However, the non-inclusion of 
a moisture-continuity equation may be more serious in regions 


of upper ridges. 


Ine inversion сазо (T.,- TA€ 0) or cooling casg. whore 
heat is lost from the atmosphere to the ocean is defined 


by an empiricism from Laevastu [12]. 
Qe = 3.0 Vio (79 - Ta) em-cal/cm"-2lıhrs (27) 


Equation (27) represents the heat flux across the 10-meter 
level. This is an empiricism requiring the 10-meter wind, 
View which will be approximated by making a frictional 
correction to the 1000 mb geostrophic wind as shown by 


equation (28): 


IK X VW Deo 








Vio = 0.5 Я. (28) 
€ 


The factor 0.5 is due to surface frictional effects and will 


be discussed in section (7), 


Before integrating equation (27) over the 0-2 layer 
the initial distribution within the layer must be determined, 
Normally cooling effects due to inversion conditions are con- 
fined to a surface layer only several hundred meters thick, 
however, since the layer mean cooling effect is desired over 
the entire 0-2 layer it will be assumed that cooling is dis- 
tributed throughout the 0-2 layer. Thus, dividing equation 


eh 


(27) by the mass of a column of air extending from the surface 
pressure Pa tc Pz, with total mass,M= (Pa, - P,)/g, gives the 
mean temperature change in a vertical column extending from 


Р, to Р, аз, 


, 


Q = Qeg (29) 
Cp СьСТА- NN 


Finally we have, ТА 


әу. РАІ AG ЊУ 


Which is the final до for the inversion case assuming 


neglible heat release from fog-droplet condensation. 


7. State-Change Parameters for the Convective Case 

Air-temperature change and wet-bulb temperature change 
equations must be derived for use in equation (26). The 
state-change equations (31) and (32) were initially developed 
by Mosby [15]. änot|1], and more recently by Boyun |l ]. 
However, in a FNWF study by Carstensen and Laevastul6 | the 
following best-fit equations were found to give hourly changes 
with a high degree of skill: 


АТА O42. (Tu - T4) - 0.10 Ме TA °C /h= (31) 


A< 


A 
t 


Equation (31) is substituted directly into the sensible- 


©15(@ш-еАл)- O.18 - Vi. Vea mb/he (32) 


heat term of equation (26) while the wet-bulb temperature 
change contributes directly to the latent-heat term. In 


25 


computing Twe it is necessary to include the small increment 
At 

of temperature found by proceeding along a mixing ratio isopleth 

from Pa, to Pa, then down a moist adiabat to Pe, (see figure 


5). Equation (33) shows this relationship. 
Да = AT Ao -XA Z 9 TE 6002 (33) 
m 


Here Дре is the one-hour height change of the condensation 
heights, determined by the height of intersection of the dry 
adiabat and the mixing ratio isopleth based upon the tempera- | 


ture - dew-point spread at the surface. 


AvZpe = Zperz ~ рг — = fa (31) 
ато A y 


Та. Та 





Figure 5. AAlug as a Function of ANE: 


Substituting equation (34) into (33) gives the local rate of 
change of the wet-bulb temperature as a function of dew-noint 


and alr-tomoornturno chaneo, ў, T 


A Tyg - АЉ Y-m (AR _ AD. (35) 
At At һ-4аЛАс AT 
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The dew-point temperature change АМ, is found by taking 
the time differential of equation (36) below | 10], 


Та-Ть 5 ош (е-е) (36) 


Here, e, is the saturated vapor pressure and e, the observed 
value, both at lO-meters, The finite-difference form of the 
time differential of dew-point resulting from equation (36) 


then follows: 


АМЬ + 0622! ел | 0.622 L_A€s (37) 
X а (Әр At дыр At 


The vapor pressure change, A& , is given by equation (32). 
Ast 
The time rate of change of saturated vapor pressure, аз 4 


is calculated using the Clausius- “lapeyron equation: 


Aes = 0,622 Le, МТА 
At ST ae (38) 


es = 6.107 exp == 5418,0 (39) 
23. lG ЖТА 


Equation (39) is an integrated form of the Clausius - Clapyron 


equation, 


There are still some unsolved parameters in the system. 
The moist adiabatic lapse rate, Jm , cannot be considered a 
constant with respect to pressure even in the limited opera- 
tional range of this study. At the condensation level, Р», 


the moist adiabat is calculated from the following [20]: 


ie? Ele 
Po + — m 


P, + ez) l e 
с. en 


All subscripts C refer to the condensation level. 





(и = (а (цо) 


Condensation parameters of temperature, Te, and pressure, 
; “Ea? 


Po, are readily calculablo from equation (11) ana (112) after 


Edson | 8]: 


Ma — NEM — (Ta- To) (41) 
(а ~ бр 


Ре = РА хо (42) 
Та 


For equation (40), the vapor pressure at the condensa- 
tion level, eç, must be calculated, Equation (43), the inte- 
grated form of the Clausius -Clapeyron equation may be used 


уч 
since air at the condensation level is saturated, and P. and 


Те аге known from (41) and (12), 


5418.0 5916, о 
=, «(6.107)exp Se apace (43) 


In like manner the vapor pressure of the sea surface, Oy? 
(from equation 32) is calpulatod since the sea surface is 
saturated (equation l/j), Also, a correction factor of 


0,98 is required for salinity effects, 


ew = (0% (ођехр | te збе] що 


Finally, the advecting wind at 10-meters, Va, must be 
derived for use in equations(31) and (32). The geostrophic 
wind at the surface is taken as equal to the 1000 mb geo- 
strophic wind. However, due to frictional interaction with 
the ocean surface, the geostrophic wind must be reduced in 
magnitude and rotated to the left, to yield the advecting 
frictional wind, Vio. For the convective case the observed 


cross-isobar angle at mid-latitudes is approximately По (10). 


Reduction of the 10-meter wind speed as expressed by 
the ratio, V/Vgo > where Vg, is the 10-пебег geostrophic wind 
speed, may be obtained from a study ру Lettau | 13], Ву 
assuming a roughness length of #„ = 0.1 cm for oceanic areas, 
the lO-meter geostrophic wind ratio for neutral conditions 
1.59 ТАРС = 0.68 lettau then allows for stability criteria, 
and for tho convoctive case a value of V/ Veo = 0.7 їз arrived 


at. Similarly, V/ Vao = 0.5 for the inversion case. 


Advection by the frictional wind is derived by first 


decomposing the vector advection as shown by equation (5): 


МАМИ ТА u DÍA 37 y la (15) 
AX д4 


Then multiplying by 0.7, and rotating the field through the 
angle K = 159 by & rotational change of coordinates leaves 


equation (l6): 


Vio Via = 0.9 (use dle Со Ж — Vago DÍA Sin 
òX 24 


— Veo DÍA Lom qo 275 Sin (46) 
> x 
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Assuming the geostrophic conditions, Uge=—— два, eL, 
equation (l6) may be written in its final form as given by 


equation (47). 


Vio’ У ТА = SH 4 cosis J (Zo т) 
+ 


- 0,9. зме” Ув 2А, 2. 
Ox ва 7 ад, (47) 


Similarly, equation (18) is derived for the advection of 


vapor pressure by tHe lO-meter wind. 


У "Мед = 07 q cosis | (2оел) 
€ 


-0). 51,15% Әғәдел | daa) (18) 
7 Ox да 24 Ба x 


8. Numerical Procedures 
All finite-difference operators used are standard cent- 


ered-differences of the type, 


V^A - MEVA J(A,B) a JUS) 


ДА «м АА = м (д. уу 


with five-point grids having a mesh distance, d = 381 km 
at 60^ latitude, Hero,m із the map factor for polar stereo- 


graph projections and d is the grid spacing. 


Ihe variables needed for input into the diagnostic omega 


equation were scaled as follows: 


30 


л 
M =M.2 
A 
Y) - Y) 2. sec”! 
а € 
GQ) =040,2 mb/sec 
A 
Wo =W,,,2* mb/sec 
D = ok cm 


^ 


T =T .2 “g 


a. 1! 
Р = P. .2 mb 
A 
е -ё .2 mb 
^ 35 
Оң -94.2 сщ“ /зес“ 
^ 
Zo 7.2 om 
: л 
9 50) ‚a gn-om*/sec? 
^ -H 
+ = 1.45812x10 singsec“' 


The following physical constants were also needed in the 


computations of CO, 


been expressed in cgs units, 


d 


1 
R 


Rv 
Ср 
L 

ўа 


38,1000, 000 
980.0 cm/sec 


0.161 x 107 
ШІ02 510! 
2800 Х 107 


cm 


2 


Шеп Xx 107 от /вес“ 
om*/sec* 


eom "/sec* 


от /вөс" 


0.9771 Х 1077 “G/cm 


fy = 0.1600 x 1077 “0/ вш 
Po = 1.213 x 10? em/cm3 
(^ = 1.055 X 10? gm/om? 
(>. = 0,919 Х 107 рпуста 


All values of those constants have 


1 


Тле scaled Helmholtz equation, which follows from 


equation (20), expressed in units of, mb - cm?/sec? 


shown on next parc, 


EN 


, 18 


PEA) Asafi d (6,0). 27 
Он 


-Зо 4 67), 422 Love а а Ое 
2 Фра" 


s48 JO) JO MN ка йг 


р 
All height values have been replaced by D-values, the 
deviation of the height field from a standard atmosphere 
height. 


The scaled stability parameter is given by equation 


(50) and is in units of cm*/sect 


N 
бн "A +: 25) 2 "4 290000 -(5- 


For computation of the stability parameter D is re- 





scaled to D= D к to prevent overflow when adding the 
standard height,Z4 —Z,&54 = 290,000 em, Stability calcu- 
lations are based on 1000-700 mb differences and throughout 
the derivation all gradients and differences at the P, and 
P, levels are assumed equal to those at 850 mb and 700 mb 


respectively. 


To prevent computational instability during the numerical 


Helmholtz solution a minimum value of stability is required. 
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This value “corresponda Tora oa ar: 
(ыа 24 < 2 = (5.-5.) 2" + 210000 | 
Equation (51) is the scaled lower boundary condition 
in units of mb/sec. 


FRE = aL 1 


am Ane 


СА ICH CAEG 


2, ^ ^ ^ 
— еда з" (51) 
2 


Here, the subscript N refers to the terrain-height dependent 
cases, 1, 2, and 3 explained in section (5). The factor uos" 


converts dynes/cm* to millibars. 


Equations for the diabatic term are computed in terms 
of heating rates, 0/0, » which are assigned to grid points 
according to the existing grid condition: convective, 


inversion, or neutral. 


Equation (52) shows the scaled heating rate for the convective 


сазе, 
A 
Са (a-i. АТА + Pe- (PP) ATi zc EU 
Cp 3600/ A€ 3600 /AT 


The factor, 3600, converts hours to second and heating rate 


• • о 
is now in units of mb -  C/sec. 


Equation (53) shows the scaled heating rate for the 


inversion case (a cooling process) also in mb - “C/sec, 


ЭЭ 


489 2 136 (05 0.341881 Du Tun)? (53) 


Ф/ — fad Є бр 


The factor (1,18) X 10) converts gm-cal/sec to dyne-cm/sec; 
K = 2h X 3600 which converts 2l-hours to seconds; 0.03 re- 
places 3.0, and permits velocity computation to be made in 


cm/sec. 


Land area grid points are masked and receive heating 
rate values of zero as do oceanic grid points which satisfy 
the neutral conditions, The remaining oceanic grid points 
receive either heating or cooling values depending on the 
prevailing condition at the point. The diabatic term is then 
computed and the scaled form in units of mb-cm*/sec” is shown 


by equation (5h). 
2. a | 
RZQ 2” (51) 
Й». | Єр 
PL 
Here, P, approximates Р, for ease of computation. 


Other equations involved in the computation of the dia- 


batic term are scaled as follows: 


(1) The air temperature change is units of °C/hr. 


Ак = 0.12 (Tio-Ta) -(0.10)2* 


(ЕР?) 
022% css Mis Ta -Sinis дү, ¿AR A] 2" 
Bad” ( ) ) у 


3h 


(2) The dew-point temperature change in units of “C/hr. 


AT ФА, абі om mm mn m 
76 Ahr ae ў 5(б.о- бо) «8 
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| а al г 3 P 
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c in 2. ОС, Tar?’ #273,16) ANE = 


The wet-bulb temperature change in units of ^C/hr. 


Ada zd аһ. бэж-д. sooo Á - 45.) 
be MM ‚(їс эб, ME At 
7: ега. . 2 
Where, (= 20 4991 Dez 4 Roa d+ 293.16) > + RL Т^ 2 4 225, Іс) 
шн л (57) 


Ср бут. 24+223. 16)“ 
(h) The condensation level temperature in units of C. 


6 ^ ^ ^ 
Te = ТА 4221 үт, т) (58) 
О. (91 
(5) ne dew-point temperature in units of до 


E ^ 
= бре es 2-64) (£9) 
» РА 


(6) The condensation level pressure in units of 


millibars, Р. мал + Per)/2 


= (60) 
лы T 
Ta - 214 292808 


A A 
Пере 16, 15 computed using values of Ta ande ‚wingönuation 56, 


N ^ 2 4 
s = РА Teg 2' +223.!6 
ША ) +27316 


zs 


^ ^ 
Here Те; is calculated using values of Таў +AT,/At and 
" 
Tp, tA T/At in equations(58), (59), (55), and (56). 


9, The Computer Program 

The Control Data Corporation 160) digital computer 
was used in this study. It has a core capacity of 32,768 
words of 8 bits each. An operational field of 1,977 
grid points forming a БІ Х 17 octagon inscribed within the 
9°N latitude circle was employed. The grid-mesh is 381 km 
true at 60° N latitude. 


Boundary conditions around the octagonal grid were 
determined by the standard subroutines used. Laplacian, 
Jacobian, and all other five-point center-difference 
operations Set the grid boundary to zero. The Helmholtz 
relaxation operation set the edge and the next interior 


border point to zero, 


Y, 
Equation (19) was solved by a two-dimensional Liebmann 
relaxation technique wherein the +1)-iterate for any point 


is givenby. 


n 
Aj - An, A| VA (AB) -Cii 
4 сар + | 


Неге, А is the over-relaxation coefficient and the residual 


аб any Step, RA, , is, 


n 
В... ХІ рај -(АВ)|- Ci 
ЯГ 2-* Bij + | 
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The over-relaxation coefficient used маз, А = 1.414, 
which allowed convergence in approximately 30 scans over 
the 1977 point grid using an initial guess field of zero. The 


| š Q) 
danverfence criterionis that the iteration ceases when € 


(пу (ин) (п) 
догіпой ру Є zA «М fT 


=; >, 


ч 1 E , NT EA (13 dl 
о аўса Ссс 


е е r -4 
Tis value corresponds to a vertical velocity of 0,3 X 10 


mb/sec for a stability of 03 X Do em”/sec” (the standard 


atmosphere Starr), 


Total computation time was approximately 1 minute and 
30 seconds with 23 seconds required for the heating term 


and 26 seconds for the boundary. condition. 


Due to the abrupt cut-off criteria used for delineating 
the three cases of diabatic heating (inversion, convective, 
neutral) the final 4%, field was smoothed with a five-point 


smoother of the form, 
- 2. 
А = A + KVA 


The smoothing coefficient k is a constant value of (1/8) 
and the field was smoothed twice. The smoothing operation 
removed small scale irregularities caused by the cut-off 


criteria of the heating term but also reduced peak values, 


А standard FIWF filtering process followed all operations 
involving the Laplacian, ‘Mis process removes all small 


scale features with wave numbers greater than 15 at latitude 
O 


Пе 


SN 


10. Results and Conclusions 
A series of four successive l2-hourly data-sots be- 

ginning with the 00Z maps of 20 April 1966 and ending with 
127, 29 April 1966 were used, Since the only diabatic 
heating mechanism introduced into this study was that aris- 
ing by "conduction" from the underlying oceanic surface, 

the results have been depicted only over the North Atlantic 
Ocean, This region has a large density of reporting ships 


and thus the reason for its selection, 


The computational omegas and their associated parameters 

for 00Z 28 April 1966 are contained in figures (6) through 
(11) which are appended. Of these, not only the diabatic 
omegas are shown, but also, the additive effect of the dia- 
batic influence as contrasted with adiabatic computations. 

In figures (12) through (17) only the final-product omegas 

of this study are shown together with the corresponding 
PNW analyzed 850 mb D-fields to serve as idontifiors, 

that is, to indicate qualitative coherency between the verti- 


cal motions and the associated motion systems. 


In the sequence of figures (6) through (11), it is of 
interest to note the pattern similarity bétween comparative 
adiabatic yertical motion computations,thosé produced here, 
and, for the same times, those, producedby FNUF, which are based 

upon the procedure described by Haltiner et a1[ 9]. The 
magnitude and position of the updraft-centers show a strong 
similarity with those of FIT/F, however, the model described 


by this paper shows larger downdraft areas extending east- 
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ward from Newfoundland, This apparent discrepancy may be 
attributed to the two differing treatments in W,,(figures 
18 апа 19), TheW,,at the continent-edge essentially becomes 


the boundary condition for the oceanic computations. 


On the other hand, however, one of the major objectives 
of this study was to test a simplified ^A, which uses only 
parameters pertaining to the standard levels of analysis 
(equation 23). In this connection, the greatest time- 
consuming aspect of the FNWF w-computation is that of 
deriving Were, and associated parameters at terrain height re- 
oulring pressure extrapolation. The general similarity in 
the fields of W¿oby the two operational procedures which 
have been discussed is evident by an inspection of figures 
(18) and (19), and this can perhaps be a justification for 


continuing the simpler Q&Locomputations. 


The diabatic effects, which were obtained for the 002 
April 28 synoptic time, depict values of А/б», So that the 
large positive center southeast of Newfoundland is realistic 
with regard to the northerly flow passing over the Gulf 
Stream. The maximum heating value for the four map periods 
over the Gulf Stream is 0/Cp = 0,067 mb - °c/sec, which 
corresponds to a layer mean heating rate of 0.3 °C/hr for the 
1000-700 mb layer, The equivalent thickness change for this 


layer is approximately 3 meters/hr. 


For comparison, Petterssen's values [27] for the 
1000-500 mb thickness change during a Similar synoptic 


Situation in late March show a maximum value of 6 to 8 


= 


meters/hr in the vicinity of the Gulf Stream. These values 
for sensible and latent heating were computed using flux 
calculations at the surface based on empiricisms similar 

to those of Laevastu, It should be noted that heating rates 
for this time period in late March are indeed greater than 


late April values. 


Finally, with regard to the diabatic term, note that 
in the area of qualitative verification no significant dia- 
batic cooling values were observed, These areas ave у те >. 
ed in extent in winter and spring, but cows d be of greater 


synoptic consequence in the summer and fall seasons. 


The effect of the diabatic term as it appears on the 
forcing function side of equation (2l) gives rise to figure 
(11), applicable at approximately 850 mb, but actually rep- 
resenting the layer-mean(Qy. As could be anticipated from 
the standpoint of heat injection into moving parcels, the 
change (inclusion of the diabatic term) has been such as 
to expand the southern rim of the trough which extends to- 
wards the southwest from Iceland. Similar aspects of coher- 
ency between the resultant w of this study and the trough 
movement and development may be traced out, feature by 
feature as one proceeds through the synoptic sequences. 
For example, the pronounced trough over the Atlantic has, 
by 002 of 29 April, become oriented North-South just east 
of Greenland, At the same time the updraft cell has taken 
on this same orientation just west of the same trough, so 


that the) field gives some confirmation of the dynamic 


lh o 


1 


processes involved in these map changes. 
e 
The model for vertical motion, as presented by this 
paper, appears to be quite sucessful and future plans involve 
inclusion of a more realistic frictional term, These more 
realistically computed omegas (extended to 500 mb) are then 
to be used for feedback to yield a prognostic z-field, hour- 


by-hour. The ultimate aim, of course, is to realize a small- 


er R.M.S. error as the motion systems progress over the ocean, 
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